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(N : ABSTRACT 
O ; 

^ ■ The dynamics of GRB jets during the afterglow phase have an important 

■ effect on the interpretation of their observations and for inferring key physical 

^ ' parameters such as their true energy and event rate. Semi-analytic models gener- 

^ ■ ally predict a fast lateral expansion, where the jet opening angle asymptotically 

^ grows exponentially with its radius. Numerical simulations, however, show a 

W ' much more modest lateral expansion, where the jet retains memory of its initial 

K ■ opening angle for a very long time, and the flow remains non-spherical until it be- 

• comes sub-relativistic, and only then gradually approaches spherical symmetry. 

^ ■ Here we suggest a new analytic model based on a new physically derived recipe 

i-{ • for the lateral expansion. We also generalize the model by relaxing the com- 

^ ■ mon approximations of ultra-relativistic motion and a narrow jet opening angle. 

We find that the new analytic model fits much better the results of numerical 

CN ■ simulations, mainly because it remains valid also in the mildly relativistic, quasi 

^ • spherical regime. This model shows that for modest initial jet half-opening an- 

^ ■ gles, 60, the outflow is not sufficiently ultra-relativistic when its Lorentz factor 

\Q • reaches F = I/60 and therefore the sideways expansion is rather slow, showing no 

Q^ ■ rapid, exponential phase. On the other hand, we find that jets with an extremely 

5 ■ narrow initial half-opening angle, of about Oq ^ 10^^'^ or so, which are still suf- 

^ ■ ficiently ultra-relativistic at F = 1/6q, do show a phase of rapid, exponential 

J> ■ lateral expansion. However, even such jets that expand sideways exponentially 

^ . are still not spherical when they become sub-relativistic. 
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Introduction 



The ultra-relativistic outflows that p ower gamma-ray bursts (GRBs) are thought to be 



colhm ated into narrow jets (for reviews see iPiran 



Piran 


2005; 


Granot 


2007; 


Granot & Ramirez-Ruiz 



201 ll ). The evidence for this is rather indirect, however, since their images are usually un- 
resolved, and in the best case (GRB 030329) the late ti r ne radio afterglow image wa s only 
marginally resolved JPrail et al.lll997l : iTavlor et al1ll997[ bood : IPihlstrom et al.lboOTh . The 
different lines of evidence for jets in GRBs include analogy to other ast rophysical rela tivis- 
tic outflow sources such as active galactic nuclei or micro-quasars (e.g. iRhoadd 119971 ) . the 
difficulty in transferring enough energy to u ltra-relativistic ejecta in a spherical explosion 
of a r nassive s tar (f or long duration GRBs; iTan. Matzner fc McKed l2001t iPerna fc Vietri 
2OO2I ; iGranotI 120071 ). extremely large isotropic eq uivalent energies in some GRBs (with 



E. 



7, ISO 



4.9MoC^ in GRB080916C lAbdo et al.l 



2OO9I ) , and an achromatic steepening o f the af- 



tergl 



o w lightcurves of some GRBs that is attributed to a jet (known as a "jet break " ; iRhoads 



1997 
1999 



199 



Sari. Piran &: Halpern 



Halpern et al 



2000 



1999 



Price et al. 



20011 ) 



Fruchter et al. 1999; Harrison et al. 1999; Kulkarni 



Therefore, there is very little direct observa- 
tional information about the jet angular structure and dynamics, which make it difficult to 
interpret GRB afterglow observations and infer from them important physical parameters 
such as the jet energy and opening angle, the external density profile, and the microphysical 
parameters of the relativistic collisionless shock powering the afterglow emission. 

Most studies of GRB jet dynamics during the afterglow phase have focused on a roughly 
uniform jet with well defined, sharp edges. We shall also focus on such a uniform jet, and 
only briefly remark on the expected relation to jets with a smoother angular struct ure (also 



1999I; 


Sari. Piran & Halpern 


I999I; 


Panaitesci 


L & MeszarosI 


1999t iKumar & Panaitescu 


2000; 


Moderski, Sikora & Bulik 


2OO0I: 


Piran 


2OO0I 


Oren, Nakar & Piran 


I2OO4; 


Granot 


2007 


) and 



numerically, using two dimensional special relativistic numerical simulations (IGranot et al. 



2001 



2010 



Cannizzo et al. 



Wygoda et al. 



2011 



2004lzhang fc MacFadvenll2009( Ivan Eerten et al.lboidlMehani fc Keppens 



van Eerten fc MacFadyenll201l[ ). as well as an intermediate ap- 



proach where the dynamical equations are integrated over the radial profile of t he thin 



shocked region, thus reducing the set of partial differential equations to one dimension (IKumar fc Granot 
2003h . 



Let us consider a uniform double-sided jet of total energy i^jet, initial half-opening angle 
^0, and initial Lorentz factor Pq. GRB observations suggest that typically ro6'o ^ 1. At 
early times, as long as P 3> 0^^^ the bulk of the jet is causally disconnected from its edge and 



(l-COS^o)-'^jet ~ 

26Q ^Ejf.i, following the spherical iBlandford fc McKed (jl976l ) self-similar solution. This early 



thus evolves as if it where part of a spherical flow with an e nergy Ei, 

3-2; ' ' 
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phase corresponds to radii R < Rj, where the jet radius Rj is defined as the radius where 
r = 1/6q for a spherical fiow with E = E[so- At R > Rj the bulk of the jet is in causal 
contact with its edge and the jet can in principal rapidly expand sideways. However, the 
degree of lateral spreading at this stage, which strongly affects the dynamics, is not well 
known. Therefore, the jet dynamics at i? > Rj are still controversial. In particular, the 
radius -Rnr at which the fiow (or jet) becomes non-relativistic, still remains uncertain. 

The Sedov length for a spherical fiow with the true jet energy, E = E-^^t (i.e. the radius 
where it sweeps up a rest mass energy equal to its own energy and becomes non-relativistic), 
Rs{Ejet), is very close to Rj. Therefore, in order for the jet to be already close to spherical 
when it becomes non-relativistic (i.e. at -Rnr), it must expand sideways very quickly and 
become close to spherical already near Rj [i.e. i?NR cannot be ^ Rj ~ i?s(-£'jet)]- This is 
indeed roughly what happens in simple analytic models, where the jet half-opening angle, 
6j, starts growing exponentially with radius near Rj, and the jet quickly becomes close to 
spherical and non-relativistic at a radius -Rim ~ (1 — ln^o)-Ri ~ (1 — In 6'o)-Rs(-E'jet), which is 
larger than Rj only by a logarithmic factoio, while at -R > -Rnr the fiow quickly approaches 
the Newtonian, spherical, self-similar Sedov- Taylor solution. 

Numerical simulations, however, suggest that most of the energy remains within the 
initial jet half-opening angle 6o until the fiow becomes mi ldly relativistic, and only then does 



the flow start to gradually approach spherical symnietrv (IGranot et al. 



2001 



ICannizzo et al. 



2004J : IZhang fc MacFadven] l2009t Ivan Eerten et al.ll2010l : iMehani fc KeppensI 120101 ). Under 
the crude approximation that the jet does not expand sideways and keeps evolving as a coni- 
cal section of a spherical fiow up until the radius where it becomes non-relativistic, the latter 
is given by -Rnr ~ -Rs(-E'iso) = 6*0 "^^^^ ^^^o- this case the fio w is still highly non-spherical at 



-Rnr, and only very gradually approaches spherical symmetry (IGranot. Ramirez- Ruiz fc Loeb 
2005h . 



This clearly shows that without lateral expansion -Rnr is significantly larger, by a factor 
of ~ ^^''^ '^V(l ~ lii^o) (which is ^ 1 for ^ 1), than if there is fast lateral expansion 
at -R > Rj. Thus, the dynamics of the fiow at small radii (-R <^ Rj) and at large radii 
(-R ^ -Rs(-Eiso)) are reasonably well known, while at intermediate radii (-Rj < -R < -Rs(-Eiso)) 
they are still controversial. For typical values of ~ 0.1 this range of radii may appear rather 
small, Rsn/Rj ~ 1 — In 6*0 ~ 3.3 for exponential lateral expansion and -RnrZ-Rj ~ Oq "^^^^ '^^ ~ 
4.6 for no lateral expansion up to -Rnr ~ -Rs(-E'iso) with k = 0. However, it corresponds to a 



^The mild discrepancy, by a logarithmic factor, between -Rnr and -Rs(-£'jct) ~ Rj likely arises from the 
fact that in simple analytic models the swept-up mass at the radius where the jet becomes spherical is smaller 
than the external rest mass within a sphere of the same radius. 
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large range in observed times (over which the corresponding afterglow emission reaches us), 
of ^obs.NR/^obsj ~ (1 — \n9f))9Q^ ~ 330 and ~ ^-(8-2fc)/(3-fe) ^ ^gg^ respectively, since the 
observed time scales as tobs ~ R/cT'^, and F decreases by a large factor (of 9q^) within this 
range of radii. 

Most simulations so far were for = 0- 2, or even wider initial jet half-open i ng an 



have p erformed simulations also for narrower initial jets, 9q 



van Eerten & MacFadven 


( 


2011) 


9a = 0.05. 0.1. 0.2. IWveoda et al. 



( 1201 if ) have found that significant lateral spreading starts when F drops below 9^ , as pre- 
dicted by analytic models, and tried to reconcile the apparent discrepancy with analytic 
models by attributing it to their small range of validity after significant lateral spread- 
ing starts (1 -C F < 9n ^ ) for the typical modest values of 6*0 used in the simulations. 



van Eerten &: MacFadvenI (120111) disagr ee with this conclusion, and we address this dispute 



in § [71 More recently, iLyutikovl ( 1201 ll ) has argued that significant lateral spreading is ex- 
pected only at a later stage, when F drops below 9q (rather than 9q as obtained in 
simple analytic models), based on an analytic consideration (which we address in § |3] and 
Appendix |X1 and find to be in error). Thus, there appears to be an ongoing debate on these 
important issues. 

Here we try to reconcile the apparent differences between the analytic and numerical 
results, in light of this recent debate. The different relevant critical radii are discussed in 
§ 121 In § [3] we discuss the recipe for lateral expansion used by analytic models, and derive 
a new recipe that takes into account the non-spherical nature of the shock driven by the jet 
into the external medium. In §|l]we construct an analytic relativistic model, which includes 
both the traditional recipe and our new recipe for the jet lateral expansion. It is also 
shown that while the region of interest and validity of the analytic model (corresponding to 
K F < ) increases as 9q decreases, 9j reaches lower values, resulting in a narrower jet at 
the time when the analytic solution becomes invalid. Because this relativistic model breaks 
down in a region of interest (both for typical GRB parameters and for comparison with 
simulations), in §|5]we generalize it so that it would be valid also at low F and high 9j, using 
two different assumptions on the accumulation of the swept-up external mass (in § 15.11 and 
§ 15. 2p . According to the results of these models (§ 15. 3p . a phase of rapid exponential lateral 
expansion exists only for sufficiently narrow initial jet half-opening angles, of approximately 
9o <^ 0.05, 0.03, 0.01 for k = 0, 1, 2, respectively. In §[6] we compare our analytic models to 
numerical simulations (with a modest = 0.2) and find reasonably good agreement (with no 
exponential lateral expansion in both cases), where the differences between the two recipes 
for the lateral spreading have a smaller effect on the agreement with numerical simulations 
compared to the generalization of the model to small F and large 9j. The implications of 
our results are discussed in § [71 
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2. The different critical radii, and two extreme assumptions for the jet 

dynamics 

Using the approximate equation for energy conservation (for T ^ 1), E ^ r^M(i?)c^, 
where M{R) is the swept-up rest mass at radius R for a spherical flow in an external density 
Pext = AR~'' (with A; < 3), and the deflnition of the jet radius Rj as the radius where T = 9q^ 
for a spherical flow of energy Eiso, we obtain 



Rj 



(3- A;)Ejet 



l/(3-fc) 



2V(3-'=)i?g(£;.^,) 



Similarly, the Sedov radius of a spherical flow with E = Eiso is 



(3 - k)Ei, 



1 i/(3-fc) 



9, 



-2/(3-fc) 







i?. = ( ^ 



2X -l/(3-fc) 




'jctj 



(2) 



Two extreme assumptions on the degree of lateral spreading, which likely bracket the 
true jet dynamics are: 1. mildly relativistic lateral expansion in the jet co-moving frame, and 
2. no lateral spreading until the jet becomes non-relativistic. Assumption 1, which is made 
in most semi-analytic models, results in exponential growth of Oj{R), until the jet becomes 
quasi-spherical and non-relativistic at 



-Rnr.i ~ (1 — ln^^o)-Ri (fast lateral spreading) 



(3) 



Assumption 2 was so far studied mainly by lGranot. Ramirez- Ruiz fc Loebl ( 120051 ) . and leads 
to 

-Rnr,2 = Rs{Eiso) = 6q'^^^^~''^ Rj (no lateral spreading) . (4) 

In this case the jet is still very far from spherical symmetry at -Rnr, and thus approaches 
spherical symmetry only after the radius grows by a factor 62, of a few or several. Moreover, 
since the radius of the Sedov- Taylor solution scales as 



RsT{E,t) ^ Rs{E) 



ct 



1 2/(5-fc) 



Rs{E) 



Ef_ 
A 



l/(5-fc) 



(5) 



and at the non-relativistic transition time, tNR,2 ~ -Rnr,2/c = Rs{Eiso) / c, the Sedov- Taylor 
radius of a spherical flow with the true jet energy is much smaller than the jet radius at 
that time, -RsT(-Ejet, ^nr,2)/-Rnr,2 ~ 6'q ^^^^ ''^ ^ 1, the flow approaches spherical symmetry 
only at the time tsph,2 when -Rsxf-E'iRt; ^sph. 2) = ^2-Rnr,,2 = fe2-Rs(-£'isn)i which cor responds to 
^sph,2/^NR,2 ~ ^62^"'^''^^ ^ 1 (see Eq. 6 of Granot. Ramirez- Ruiz fc Loeb 2005 ). Note that 



this is much smaller than the factor (62) by which the radius grows over the same time. For 
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log(t) 

Fig. 1. — A schematic figure showing the evolution of the jet radius R = R\\ (i.e. its 
extent along its symmetry axis) and lateral size R± as a function of the lab frame time t 
for two extreme assumptions on its degree of lateral spreading: (1) mildly relativistic lateral 
expansion in the jet co-moving frame, and (2) no lateral spreading until the jet becomes 
non-relativistic. The jet becomes spherical when R\\ and -R_l become equal, which occurs 
well after the jet becomes non-relativistic, and then joins the Sedov- Taylor solution. 

assumption 1 similar arguments imply tsph.i/^NR,! ~ (1 where bi < 62 

can be expected. 

Fig. [1] shows the jet radius R = R\\ (i.e. its extent along its symmetry axis) and lateral 
size R± as a function of the lab frame time t for these two assumptions. The region where 
the dynamics for these two extreme assumptions differ is basically where the dynamics are 
most uncertain, and corresponds to the range of radii i?s(-^jet) < R < 62-Rs(-^iso) (i-e. a 
factor of fn ~ 6'q'^''''^~'^^62 in radius), and (lab frame) times i?s(-^jet)/c < t < tsph2 ~ 
9,'bt'^/'Rs{E,,o)/c (or a factor of /, ~ ff^^' ~ ^-(5-.)/{3-.)^(5-.)/2 
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Altogether, the ordering of the different radii is 

Rs{Ejet) ~ Rj < -RnR,1 < -Rs(-E'iso) = -RnR,2 , (6) 

or 

-Rs(-Ejct) q _ In g^)"^ "^^^'^ ~ g-2/(3-fc) -Rs(-£^iso) ^ g-2/(3-fc) -RnR,2 



p . ^ p . ^ p . ^ p . 

J J J J 



3. Analytic recipe for lateral expansion 

The "traditional" basic underlying model assumpti ons used for the analytic modeling of 
relati vistic jet dynamics during the afterglow phase (e.g., Rhoads|[l999 : Sari. Piran fc Halpern 



19991 ) are (i) a uniform jet within a finite half-opening angle 6j with an initial value Oq that 
has sharp edges, (ii) the shock front is part of a sphere at any given lab frame time t, (iii) 
the outer edge of the jet is expanding sideways mildly relativistically, with ~ 1 in the 
local rest frame of the jet (where quantities are denoted with a prime), (iv) the jet velocity 
is always in the radial direction and 6*^ -C 1. Under these assumptions, the jet dynamics are 
obtained by solving the ID ordinary differential equations for the conservation of energy and 
particle number^! 

The lateral expansion speed in the lab frame (i.e. the rest frame of the central source and 
the external medium) is (3g = uq/T = u'g/T, where ug = Tl3g is its lateral component of the 4- 
velocity (which is Lorentz invariant, so that u'g = ug), while Ur = is its radial component. 
Primed quantities are measured in a frame moving at i3r f in the radia l direction, so that (3'^ 



and (3' = [1 — (F') ^]^/^ = (3g. The usual assumption (lRhoadslll999l : ISari. Piran fc Halpern 
19991 ) is that ~ 1, which corresponds to 

Pe-^. (8) 

As is shown in the next section, I3g ~ dOj/dlnR directly determines the jet lateral expansion 
rate in the lab frame. 

Here we derive a new physically motivated recipe. It relies on the fact that for any 
shock front, with an arbitrary shape, the local velocity vector of the material just behind 



^For the adiabatic energy conserving evolution considered here, the equation for momentum conservation 
is trivial in spherical geometry, and does not constrain the dynamics. For a narrow {9j <^ 1) highly relativistic 
(r ^ 1) jet, the equation for the conservation of linear momentum in the direction of the jet symmetry axis is 
almost identical to the energy conservation equation. When the jet becomes sub-relativistic the conservation 
of energy and linear momentum force it to approach spherical symmetry, and once it becomes quasi-spherical 
then again the momentum conservation equation becomes irrelevant. 
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the shock front as measured in the rest frame of the upstream fluid ahead of the shock (i.e. 
the lab frame in our case), 0, is normal to the shock front (i.e. in the direction of the shock 
normal, h, at that location; iKumar fc GranotI l2003l ) . namely 



P = n . 



(9) 



A simple way to understand this result is that as each fluid element passes through the shock 
it samples only the local conditions (and is not aware of the large scale or global shock front 
geometry) and locally the shock normal is the only preferred direction in the upstream rest 
frame (e.g., the pressure gradients that accelerate the fluid element are in the — n direction 
and thus accelerate it in the n direction). For an axisymmetric shock (with no dependence 
on the azimuthal angle 0), Eq. immediately implies that the angle a between the shock 
normal, n, and the radial direction, f, which is defined by cos a = n ■ f = /3 ■ f, satisfies 



tan a 



1 OR 
R~de 



d\nR 

de 



(10) 



where 9 is the polar angle measured from the jet symmetry axis. Since R ~ /3ct, we have 
dhiR/dO ~ dlnP/dO = r~'^dlnu/d6 ~ —l/T^AO, where is the angular scale over which 
u varies significantly, and we have assumed that u decreases with 6, as is usually expected. 
Since for F ^ 1 and a ^ 1 we also have (3^ ~ 1, Eq. (fTOj) implies that ~ 1/F^A^. For a 
roughly uniform jet of half-opening angle 6j we have A6 



6j, and therefore 



/3e 



1 



1 



'IV. 



which is our new recipe for lateral expansion. 



Eq. (ITT]) was fi rst derived in the context of GRBs bv lKumar &: GranotI (120031) . Recently 
i t was rederived by iLyutikovl ( 2011 ). based on an earlier work by IShapirol dlOToj ). ILyutikovl 
(I2OIII ) has argued that Eq. (ITT]) implies a negligible lateral expansion as long as F > 
1 / suggesting that with this model one obtains a slow sideways expansion, as seen in 
the numerical simulations. However, as we show later, this formula results in a slower lateral 
expansion (compared to the usual recipe, i.e. Eq. [8]) only as long F > l/9j (the standard 
condition for the onset of significant lateral expansion) , but once T < l/6j this formula leads 
to a faster sideways expansion. We also show later that other factors, namely the break 
down of the ultra-relativistic and small angle approximations, are the main cause for the 
discrepancy between the existing simple analytic models and the numerical simulations. For 
completeness we discuss the details of Lyutikov's and Shapiro's work in Appendix A. 



- 9 - 



4. A simple relativistic model 



We turn now to compare the traditional recipe for the lateral expansion speed, /^g ~ l/F 
(Eq. [8]), with our own new simple recipe, (3$ ~ l/r^6'j (Eq. [H]), which was derived in the 
previous sectio n. These recip es are implemented here within the semi-analytic model for the 
jet dynamics of lGranotI (120071 ). The main results are provided here and we refer the reader to 
that work for more details on that model. Broadl y similar serni-analytic models, with some 
varia t ions, were used earlier by other authors (e.g.. 



1999 



2000 



..^^ ^ „ 

Panaitescu & Meszaros 


1999 


Oren, Nakar & Piran 


2004 


)• 



Rhoads 



1997 



3 



1999I: ISari. Piran fc Halperru 



2OOOI : 



Moderski. Sikora fc Bulik 



The lateral size of the jet, R±, and its radius, R = R\i, are related by R± ~ 6jR. The 
evolution of R± is governed by 

dRi_ ^ OjdR + (3ecdt ^ {9j + (3e) dR , (12) 

and therefore 



ddi 



1 



W = n) 



dlnR 



(30 



(13) 



(«^~1), 

where we have conveniently introduced the parameter a that enables us to analyze these two 
different recipes together. 



The external density is assumed to be a power law in radiuqj, pext = AR-''. The total 
swept-up (rest) mass, M{R), is accumulated as 



dM 
~dR 



27r(^,/?)Vext(i?) = 2nAR''-'e^AR) 



(14) 



where the factor of 2 is since a double sided jet is assumed. As long as the jet is relativistic, 
energy conservation takes the form Ejet ~ F^Mc^, which implies that Md^T"^) = —T'^dM, 
and 

c^r T dM . „,_,„o.„. T{R) 



dR 2M dR ~^^^'"'^?^^)M(i?) 



(15) 



One can numerically integrate equations (fT3|) . ( Il4|) . and ffTSj) thus obtaining 9j{R), M{R), 
and T{R). Alternatively, one can use the relation E^et ~ F^Mc^ (energy conservation) which 



■^We consider here and throughout this review only fc < 3 for which the shock Lorentz factor decreases with 
radius for a spherical adiabatic blast wave during the self-similar stage of its evolution ( Blandford fc McKeel 



19761) . 
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reduces the number of free variable to two, and solve equations f[T3|) and f[T5|) . Changing to 
normalized dimensionless variables 9 = 9j/9o, 7 = T9q and r = [(3 — k) /2Y^^^~''^ R/ Rj, gives 

H9 

^ = r-^-f-'-''(r)9-\r) , (16) 
dr 

^ = -r^-N^(r)^^^(r), (17) 
dr 

where the initial conditions at some small radius tq ^ 1 (just after the deceleration radius) 
are 

0{ro) = 1 , 7(ro) = sj^r,^'-'^^' . (18) 
Note that by definition, 7^ = r9j. Eqs. f lTC]) and f[T7|) imply 



(ir r(76')" r(76')" 



For r <^ 1 the second term on the r.h.s of Eq. flTOj) dominates, implying (7^)^ ~ ^^r'^"^, 
which is consistent with Eq. f lTSj) . This suggests that the two terms become comparable at 
r ~ Tc that is given by 

3 _ i;\ (3+a)/[(l+a)(3-fc)] 

— ) ■ 

While Tc > 1 for A; < 1, it can reach very low values (r^ ^ 1) as A; approaches 3. We 
are interested here mainly in > 2, for which Tc ~ 1 still approximately holds. We do 
note, however, that the lower values of Tc for higher values of k result in an earlier onset of 
significant lateral expansion for such higher /c-values. Now let us examine what happens at 
r ^ Tc ~ 1. If we assume that the first term becomes dominant then Eq. f[T^ would imply 
7^ ~ [(1 + a) lnr]^/(^+"\ which in turn implies that the second term would be dominant 
(since (^^)3+a^3-fe ^ ^^-^ _|_ i]2 7-](3+«)/{i+«)7'3-'= ^ 1)^ rendering the original assumption 

inconsistent. The same applies if the opposite assumption is made, that the second term 



is dominant (in this case 'y9 ^ y ^_k^{k 3)/2 -^j-^j^j^ implies that the first term would be 
dominant, (-y6l)3+a,^3-fc ^ ^3^j ^(A;-3)(i+a)/2 ^ -^y rjj_^jg ijj^piigg f;[^gx the two terms 

must remain comparable, implying 7^ ~ j,(fc-3)/(3+a)^ similar conclusion can be reach by 
taking the ratio of equations ( IT6l) and (|T71) which implies that 

d{9'^+^) = r''-^d{-f-'^-^) . (21) 

A more careful examination shows that they must cancel each other to leading order, and 
the first two leading terms for r ^ 1 are given by 

76* ^ ^(fc-3)/(3+a) I 3 ~ ^ ^(fc-3)(2+a)/(3+a) _ ^22) 

(3 -|- Q.)^ 



Substituting equation (122|) into equations f|T6|) and f|T7|) yields 



and 



ciln^ 
dlnr 
din 7 
(ilnr 



,(3-fc)(l+a)/(3+a) _ (3 - A:)(l + o) 

(3 + a)2 

_^{3-A;)(l+a)/(3+a) _ ^(3 - k) 



(3 + a) 



(3-fc)(l + a) 

^ 6 r (3+")^ exp 



]^ 2(3-fc) 





(3 + a) 



(3- A;)(l + a; 



(3-fc)(l + a) 
7- (3 + a) 



(3 + a) 



(3-fc)(l + a) 
f {3 + a) 



{3-k){l + a) 

where the normalization coefficient b is determined numerically. For r ^ 1 we have 

'3 - 



(23) 
(24) 

(25) 
(26) 



'j9 



3 — k 

r 2 + 



-2-a a(3-fc) 



(3 + a) 



(27) 



Fig. [2] shows the results of our model in terms of the normalized jet half-opening angle 
6 = Oj/60 and Lorentz factor 7 = TOq (as well as their product, = T6j) as a function of 
the normalized radius r = [(3 — k) /2]^^^^~''^ R/ Rj . The results are shown both for a uniform 
external medium [k = 0), which is the main focus of this work, as well as for a stellar wind 
{k = 2; this is included mainly for completeness and is only briefly discussed in § [7]). The 
dynamical range in this figure is unrealistically large, and it is shown mainly in order to 
demonstrate the properties of this solution, and show how well our analytic approximation 
for r > 1 works (the dashed green lines in the middle and bottom panels, which are practically 
on top of the numerical results). The excellent agreement between our semi-analytic results 
(the numerical solution of Eqs. [16] and [17]) and analytic formulas (Eqs. [25] and [26]) shows 
that our analytic results (including Eqs. [22] and [27]) can be safely used in order to analyze 
the result of this model. This good agreement was also used in order to find the exact 
values of the numerical coefficient b that determines the normalization for 6 and 7., which 
were found to be b{k = 0,a = 0) ^ b{k = 0,a = 1) ^ 0.60, b{k = 2,a = 0) ^ 0.395 and 
b{k = 2,a = 1) ~ 0.45. Our new recipe for the lateral expansion speed (Eq. [11]) results in 
a slower initial lateral expansion compared to the old recipe at r ^ 1, where TOj = ^6 ^ 1. 
However, at larger radii, r > 1, where F^^ < 1 it results in a faster lateral expansion. 

Figure [3] shows similar results for a uniform external medium {k = 0) and for three 
different values of the initial jet half-opening angle, 6q = 0.05, 0.1, 0.2. Since the dynamical 
equations (Eqs. [16] and [17]) involve only the normalized variables 6, 7 and r, and the 
initial conditions (Eq. [18]) for 6 and 7 depend only on the initial normalized radius tq, the 
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Fig. 2. — The jet dynamics according to our relativistic analytic model (see text in § H] 
for details), for either a uniform {k = 0) or a wind-like, stratified {k = 2) external density 
profile, and for either the old (a = 0) or our new (a = 1) recipe for the jet lateral expansion 
speed. The dynamical range in this figure is unrealistically large, and it is shown mainly 
in order to demonstrate the properties of this solution, and show how well our analytic 
approximation for r > 1 works (the dashed green lines in the middle and bottom panels, 
which are practically on top of the numerical results). 
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r - R / Rs(Eje,) 

Fig. 3. — Similar to Fig. [2] but only for a uniform external density {k = 0) and for three 
different values of the initial jet half-opening angle: 6^0 = 0.05 (green), 9q = 0.1 (red) and 
6*0 = 0.2 (blue). The old (a = 0) and new (a = 1) recipes for the jet lateral expansion are 
sown by dashed and solid lines, respectively. In the top two panels the lines for different 
and the same a values coincide (see text for details). The values of Rs{Eiso) are indicated in 
the bottom panel for reference. 
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lines for these normalized variables in the top two panels for the different Oq values exactly 
coincide0 The two bottom panels show the un-normalized quantities 6j and T for our three 
values of 6*0- In the bottom panel we have added for comparison the Sedov radius, i?s(i?iso), 
for a spherical flow with the same isotropic equivalent energy the jet started with. We define 
i?NR for our model as the radius where formally F = 1 (at which point this model clearly 
breaks down). Figure H] is similar to Fig. |3]but the jet radius R is normalized by the radius 



-2/(3-fc) 







Rj instead of [(3 - k)/2]-^/^^-''^Rj = [(3 - fc)/4]-V(3-fc)i?g(^. 



'jetj 



Fig. shows -RnrZ-Rj and 6'j(_Rnr) as a function of 6*0. It can be seen that -Rnr de- 
pends on 9o only logarithmically (as can also be seen from Eq. [26]), while Rs{Eiso) / Rj = 



9, 



-2/(3-fc) 







is simply a power of Oq. It is also evident that Oj{R 



NRj < 1 for 9o -C 1, and 
its value increases with 6*0 (while 0j{Ri^^)/9o decreases with 6*0). This can also be seen 
from Eq. fl221) . using the leading order term in r and the definition r(rNR) = 1, which im- 



ply that %(rNR) 



'(3-fc)/(3+a) 



NR 



while tnr (or -Rnr) decreases (logarithmically) with Oq. 



For k 
of r = 

the jet also starts to spread sideways significantly at smaller values of r, of the order of 

f3+a)/[(l+a)(3-A:)] 



= 2 the jet becomes non-relativistic and the model breaks down at smaller values 
[(3 — k) /2Y^^^~^'> R/ Rj compared to = 0, which is consistent with the fact that 



(¥)' 



However, we are primarily interested here in /c = 0. 



The model breaks down when F drops to 1 (or even slightly earlier). As can be seen 
from Figs. [2HU it breaks down earlier for larger 9q values, and its region of validity (especially 
at i? > Rj) decreases as 6q increases. In particu lar, for the value o f 6*0 = 0.2, which was 



most widely used so far in ri umerical simulations (IGranot et al.l I2OOII : IZhang fc MacFadyen 



2OO9I : Ivan Eerten et al.ll2010l while an even larger value of 9^ 



20= 



0.35 rad was used in 



some works - iMeliani &: Keppensll2010l : Ivan Eerten et al.ll201ll ). this dynamical range is very 



narrow, and the asymptotic exponential growth of 9j with R is not reached before the model 
breaks down (at F < 1.5 — 2 or 9j ^ 0.5 — I)- Even for 9o = 0.05, which w; as used in the most 



recent simulations (IWygoda et al.l I2OIII : Ivan Eerten fc MacFadyenI 1201 ll ) and is at the low 
end of the values inferred from afterglow observations, the asymptotic exponential regime 
is only barely reached before the model breaks down (in agreement with the conclusions 
of IWygoda et al.l 1201 ll ). Note that in this (limited) region of validity of this semi-analytic 
model our new recipe might still result in smaller or comparable values of 9{r) (i.e. of 9j for a 
fixed 6*0, at a given radius for a fixed E^et) compared to the old recipe. The discussion about 
when this model breaks down is expanded in § [6l where we compare the analytic models 
to numerical simulations. Because of this important limitation of our relativistic analytic 
model, in the next section we generalize it so that it would not break down when the jet 



^This is since the same value of tq — 0.4 was used, but in the limit ^ 1 the dependence of the solution 
on ro goes away at r ^ rp. 
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R / Rs(Eiso) 



Fig. 4. — Similar to Fig. |3] but shown as a function of the jet radius R normahzed by 

i?s(i?iso) = R^e-^'^'-''^ instead of of [(3 - k) /2]-^l^^-^) R, = [(3 - fc)/4]-V(3-fc)i?s(i^.^J, 
where Eiso ~ Ejct^/^o is the isotropic equivalent energy in the jet, while E-^ct is its true 
energy. 
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Fig. 5. — Upper panel: the non-relativistic transition radius for our analytic relativistic 
model, i?NR, defined by r(_RNR) = 1, normalized by Rj, as a function of Oq. For comparison, 
i?s(-Eiso) is also shown; the two radii, i?NR and i?s(-Eiso), become similar at ~ 1 but are 
very different for 1- Lower panel: the value of the jet half-opening angle, 6j, at -Rnr? 
where our simple analytic relativistic model breaks down. 
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becomes sub-relativistic or wide. 



5. Generalized models valid for arbitrary T and 6j 

In order to avoid the breakdown of the model at small Lorentz factors F or large jet 
half-opening angles 6j, we construct here simple generalizations of the analytic model studied 
in the previous section, which do not require the jet to be very narrow {6j ^ 1) or highly 
relativistic (m ~ F ^ 1). Two variants are introduced, named the trumpet model (in § 15.11) 
and the conical model (in § 15.21) . according to the shape of the region from which the external 
medium is assumed to have been swept up by the jet (before it becomes spherical). 

The rate at which the jet half-opening angle, 6j, increases depends on the lateral velocity 
at the edge of the jet, /Sg, as d9j = f3gcdt/R = {f3g/ f3r)dR/R, or 



dOj _ (3e 



dlnR (3r 



(28) 



A crude approximation for for the comoving 4- velocity of the lateral expansion {u'g), which 
would roughly correspond to the sound speed both in the relativistic and in the Newtonian 
regimes, is ~ /3 = ^(1 + -u^)"^/^. This would modify the traditional recipe to (3g = 
u'g/T ~ (3/T = u/{l + u^) or /3g//3r ~ /3e/(3 ~ 1/F = (1 + u'^)-^/\ In our recipfl /3e//3r = 
-d\nR/de ~ -d\nl3/de ~ -T^^dlnu/dO ~ l/T^AO ~ l/T^ej = + u^)9j]. Therefore, 
just as before, we still have 



de, 



J _ 



dlnR f3r T^+^e^ 



= n), 

(29) 

K ~ 1) • 



5.1. The "trumpet model" 

In this model we follow the usual assumption that the external rest mass is swept-up 
by a working area consisting of the part of an expanding sphere of radius R within a half- 
opening angle Oj{R). Thus, the total swept-up (rest) mass, M{R), for a double-sided jet is 



^Note that we use d\au/dO = u^^du/dO ^ —1/Ad since the 4-velocity u, unhke (3 or T, generally varies 
significantly with 6 both in the relativistic and in the Newtonian regimes, so that du/dO ^ —u/A9 in both 
regimes, while dT /dO ^ —T//S.9 only in the relativistic regime and 0(3/89 ~ —/3/A9 only in the Newtonian 
regime. 
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accumulated as 

^ ^ [1 - cos^j (i?)]47ri?Vext(i?) = [1 - cos%(i?)]47rAi?2-fc . (30) 
dR 

Energy conservation takes the approximate form E^et ~ u'^Mc^, implying Md{u'^) = —v?dM, 
and 

du u dM 2tiA(? ^ „ , s., -j., s. , ^ 

IR = -mlR = -^R^-'\l-cosOAR)]uHR) . (31) 

Thus, in terms of r = [(3 — k)/2Y^'^^^''^R/Rj we have 

- JTTW^ . J - -^^-'-'('O^ll - cos MOI . (32) 

where the initial conditions at some small radius Rq ^ Rs{Ejct) ~ Rj (just after the decel- 
eration radius), corresponding to tq, are given by 



5.2. The "conical model" 

Here we note that the usual assumption that leads to Eq. (l30i) neglects the external 
matter at the sides of the jet. Because of this, when eventually 9j reaches 7i/2 at -Rgph and 
is thus assumed to be fully spherical, the amount of swept-up external rest mass at -Rsph 
calculated according to Eq. (!30l) will be significantly smaller than that originally within a 
sphere of the same radius. Therefore, here in the conical model we adopt an alternative 
approach of using for the rest mass of the swept-up matter, that originally within a cone of 
half-opening angle dj, 

4tt 

M{R) ^ [1 - cos % {R)] ttAR^ ■ (34) 

(3 — k) 

This still has the drawback of assigning the same Lorentz factor to all of the swept-up 
external matter, even though that at the sides of the jet should have a significantly smaller 
4-velocity than that near the head of the jet. Using a slightly different normalized radius, 
rg = R/Rs{Ej^t) = 2^/^^-^^RlRj = [4/(3 - k)Y^^^~''^r, energy conservation (Ejet ~ m^Mc^) 
and Eq. (!29|) imply 

r"^^-''^/^ dO- 1 

Vl-cos%(rs) ' c^lnrs [1 + r^-=^(l - cos0,)-i](i+-)/2^« ' ^ ^ 
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where the initial conditions at some small radius Rq <^ -RNR,sph(-E') ~ Rj, corresponding to 
rs,o, are given by 



{3-k)/2 



Ojirsfl) = 00 , M(rs,o) = ' „ • (36) 

a/I — COSC^o 



5.3. Results for the generalized models 

Figures El [7] and M depict a comparison of these two models with the relativistic model. 
All three models agree at early times, while the jet is still highly relativistic, narrow and 
hardly expanded sideways. The approximations of our relativistic model hold well at this 
stage and the difference in the swept-up mass between the trumpet and conical models is 
still very small. At later times, however, the three models show a different behavior. The 
main effect of the relaxation of the small 6 and ultra-relativistic approximations is that for 
typical values of ^ 0.05 the region of exponential growth of 9j with R largely disappears, 
and is replaced by a much slower, quasi-logarithmic growth. This can most clearly be seen 
by comparing the results of the relativistic model (from § HJ solid lines in Figs. |6] and |8l 
and green, red or blue lines in Fig. [7]) and the trumpet model (from § 15.11 dot-dashed lines 
in Figs. OandlHl and black, magenta or cyan lines in Fig. [7j). These two models share the 
same assumption on the accumulation of the swept-up external medium, and differ only by 
relaxing in the trumpet model the requirements of F 1 and 6j <^ 1. The results of these 
two models are very close at early times while F ^ 1, but diverge as F becomes more modest 
and the simple relativistic model reaches the exponential regime. This can also be seen in 
Fig. ISlthrough the fact that 9j{r) for the two models start diverging when 9j becomes modest 
and the small angle approximation breaks down. 

The main difference between the trumpet and conical models is that for the conical 
model the swept up mass at a given radius R is larger than for the trumpet model, resulting 
in a smaller F and therefore also a larger 6j, i.e. a faster evolution of 6j and F with R. Since 
the larger swept-up mass comes from the sides of the jet, it becomes important only once 
the jet starts expanding sideways significantly, which occurs at (see Eq. |20]). This can 
be clearly seen in Fig. [HI where the dot-dashed (or solid, which practically coincide at early 
times) and dashed lines, for the trumpet (or relativistic) and conical models, respectively, 
start diverging near r^. Note that this remains valid for all fc-values, while decreases 
with k. Fig. [8] also shows that for sufficiently small values of 9o, roughly 9q ^ 0.05 for 
k = and even somewhat smaller 6q values for larger k values, there is still a phase of 
quasi-exponential lateral expansion for < r < r{6j ~ 10^*^'^) or 1.5^o ^ dj ^ lO""'^. For 
such extremely small values of the difference between the conical and trumpet models 
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Fig. 6. — Comparison between our relativistic [solid lines), trumpet [dot- dashed lines) and 
conical [dashed lines) models, where all models use our new recipe for the lateral spreading 
of the jet [a = 1), and for a uniform external medium [k = 0). Results are shown for three 
different values of the jet initial half-opening angle: Oq = 0.05 (in green), 6'o = 0.1 (in red), 
and 6q = 0.2 (in blue). For reference we also indicate the values of TOj = 1 in the top panel, 
some relevant values of 9j in the two middle panels, as well as the values of i?s(-E'iso) and 
r = 1, 2 in the bottom panel. 
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R / Rs(Eiso) 

Fig. 7. — Similar to Fig. [6] but shown (1) only for our relativistic model (green, red, and 
blue lines for 9q = 0.05, 0.1, and 0.2, respectively) and trumpet model {black, magenta, and 
cyan lines for 9o = 0.05, 0.1, and 0.2, respectively), (2) for both the old recipe (a = 0; 
dashed lines) and our new recipe (a = 1; solid lines) for the jet lateral expansion, and 
(3) as a function of the jet radius R normalized by i?s(-E'iso) = RjO^'^^^^'''^ instead of of 
[(3 - fc)/2]-V(3-^0/?^. = [(3 _ k)/A]-y('->'^RsiE^et). 
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r » R / Rs(Eje,) 

Fig. 8. — Comparison between our relativistic {solid lines), trumpet {dot-dashed lines) and 
conical {dashed lines) models in terms of the evolution of the jet half-opening agle dj with 
the normalized radius r, for = 0, 1, 2 {top to bottom panels), where all models use our 
new recipe for the lateral spreading of the jet (a = 1). Results are shown for logio(6'o) = 
—3, —2.5, ... , —0.5 {using different colors) while the values of 6o = 1, 7r/2 and the critical 
radius Tc (given by Eq. [20], where the lateral spreading is expected to become significant) 
are shown for reference. 
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t/tj 

Fig. 9. — The jet dynamics according to our different analytic models, for = 0.1 and 
k = 0. We show the jet 4- velocity, u {upper panel) , half-opening angle, 6j {middle panel), as 
well as its normalized radius ry = r and lateral size = r sin 6j {bottom panel), as a function 
of the normalized lab frame time, t/tj, for our relativistic {green lines; until it breaks down 
at r ^ 1), trumpet {red lines) and conical {blue lines) models. The solid and dashed lines 
are, respectively, for our new recipe (a = 1; Eq. [TT]) and the old recipe (a = 0; Eq. [S]) for 
the jet lateral expansion. 
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10"^ 10"^ 10"^ 10"^ 10"^ 10"^ 10° 10^ 10^ 

Tlos / "'"losj 



Fig. 10. — Similar to Fig. |9]but as a function of the observed time, Tios, at which photons 
from the front of the jet reach an observer located along its symmetry axis, normalized by 
its value at the jet break time Tiosj. 
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becomes large during the exponential sideways expansion phase, where the lateral expansion 
is faster in the conical model. We note, however, that such extremely narrow initial jet 
half-opening angles are below the smallest values that have so far been reliably inferred from 
GRB afterglow modeling, so that they might not be very relevant in practice. 

Figs. [9] and [10] show the jet dynamics according to our different analytic models, for 
^0 = 0.1. It can be seen the the differences between the various models are rather small 
until the point where our relativistic model breaks down. The behaviour of the jet radius 
{R = -R||) and lateral size (-R_l) as a function of the lab frame time (t) shows a lot of similarities 
to the analytic expectations (compare the bottom panel of Fig. |9] to Fig. [1]). This, again, 
demonstrates that our new recipe for the lateral spreading of the jet results in slower lateral 
expansion compared to the old recipe (and is closer to assumption 2 of no lateral spreading 
- dashed red lines in Fig. [T]) at early times when F > 6j but faster lateral expansion at late 
times when F < 6j (i.e. closer to assumption 1 of fast lateral spreading - solid blue lines in 
Fig. [I]). 



6. Comparison with numerical simulations 

We turn now to a comparison of our analytic models with the results of full 2D special 
relativistic hydrodynamic simulations. To do so one needs first to define which quantities 
should be compared. This, however, is not unique and can be done in different ways. For 
the 4- velocity, u, and as one (out of a few) reference value for the jet half-opening angle, 6j, 
we use the weighted mean over the energy E in the lab frame (excluding rest energy) of u 
and 9, respectively, 

, , f dEu f dE9 , , 

For the jet radius (or parallel size, = K) and lateral size (-R±) we use: 

, ^ > , > f dE z ,, ,. 2,, 2 \ dE rcvi , x 

= 1,z)e = ^jj^ , {Ri.) = {x)e = {y)E = - (rcyi)E = - J . (38) 

These averages reduce to R\\ = R± (or {R\\) = {R±)) for a spherical flow. 

In order to perform a proper comparison to our analytic models, we need to calculate 
similar averages for our jet, which at any given time is the part of a thin spherical shell 
within a cone of half-opening angle 6j. Thus, the radial integration drops out and we are 
left only with an integral over between fij = cos9j and 1, 

R\\ _ £, dfxfx _ sin^ R± _2 dfi \fl^ ^ 26 ^ - sin(20j) 



R J'^^dfi 2(1 - cos 0,)' R 71 J^^dfi 27r(l-cos^^,) 



(39) 
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We can see that R\\ = R± for 6j = 7r/2, as it should. 

Similarly, one can calculate {9) e a.s a. proxy for 9j in our models, 

I^^ arccos(/i) j^^ dO 9 sin 9 sin 9j - 9j cos 9j . , 

^^^^ -JlJ^ 1-COS0, - 1-COS0, • 

This shows that (^)e ^ (2/3)^^ for < 1, while = 1 for = 7r/2 (which is the 
value for any spherical flow, also one with a radial profile) and (2/3)6'j < {9)e < {2/11)9 j for 
< 9j < 71/2. One can also calculate the angle out to which a fraction / of the energy is 
contained (or the energy 100/ percentile), 

6^/ = arccos [1 - /(I - cos6'j)] , (41) 

and compare it to the corresponding value from the numerical simulations. 

Figures [TT] and [T2l show a comparison (for k = and ^0 = 0.2) between t he results of our 



analyt ic models and of 2D special relativistic hydrodynamic simulations (from lDe CoUe et al. 



2011al lb[). when quantifying all of them as discussed above. As can be seen from Fig. [TTl our 



models provide a reasonable overall description of the full hydrodynamic simulations, and 
thus appear to catch the basic underlying physics, despite their obvious simplicity. 

Fig. [12] shows three different ways of quantifying the jet half-opening angle, namely the 
weighted mean over the energy, {9)e {bottom panel), and two different energy percentiles, 
9o.75 {middle panel) and 6'o.95 {top panel), i.e. the values of 9 up to which 75% and 95% of the 
energy, respectively, is contained. It can be seen that 6*0.95 provides the best match between 
our analytic model and the numerical simulations. For 6^0.75 or (^) e the match is not as good 
(though even then the difference is not very large) . This might be attributed to the fact that 
our analytic models assume a uniform energy per solid angle, e = dE/dQ, within the jet 
opening angle {9 < 9j), while in practice (or in the numerical simulations) it drops towards 
the outer edge of the jet. The drop in e from the jet axis towards its edge causes both smaller 
values of {9)e and smaller values of 9f for the lower energy percentiles (or /-values) relative 
to a uniform jet with the same 9f for a large energy percentile (or /-value; e.g., / = 0.95 in 
our case). The results for our new recipe for the jet sideways expansion are somewhat closer 
to the numerical simulations compared to the usual recipe for {9)e and 6*0.75, while the usual 
recipe is perhaps slightly closer for 6*0.95. 

Both the analytic models and the numerical simulations show that the flow becomes 
spherical more than a decade in time after it becomes sub-relativistic (which may be quan- 
tified as the time when {u)e = 1)- This can be attributed to the fact that once the flow 
becomes sub-relativistic its sound speed quickly drops, and so does the rate of lateral expan- 
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Fig. 11. — Comparison, for 6*0 = 0.2 and A; = 0, between our analytic mod els (thin lines) 
and th e results of 2D special relativistic hydrodynamic simulations (from iDe CoUe et al. 
2011al lbl) o f a jet with initial conditio ns of a conical wedge of half-opening angle 6o taken 
out of the iBlandford fc McKed (119761 ) self-similar solution (thick dot-dashed black line), in 
terms of the jet 4- velocity (u), half-opening angle {6j) as well as normalized parallel (ry) and 
perpendicular (r_|_) sizes. The green, red and blue lines are for our relativistic, trumpet, and 
conical models, respectively. Thin solid lines are for our new recipe for lateral expansion 
(a = 1) while thin dashed lines are for the old recipe (a = 0). 
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Fig. 12. — Similar to Fig. [TT]but for three different ways of quantifying tfie jet half-opening 
angle, 9j. The top panel and middle panel show two different energy percentiles, ^0.95 and 
^0.75, respectively, i.e. the values of 9 up to which 95% and 75% of the energy is contained. 
The bottom panel shows the weighted mean over the energy, (6) e- 
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sion. Moreover, as the flow gradually becomes more spherical the lateral gradients become 
smaller, which makes the flow approach spherical symmetry more slowly. 

The numerical simulations show that 6f corresponding to lower energy percentiles (or 
/-values) approach their asymptotic values for a spherical flow at later times. This shows 
that the transfer of energy to larger ^-values is the slowest near the center of the jet and 
larger near its edges, which may in turn be attributed to the lateral gradient (say, of e) in 
the jet, which are smallest near its center and largest near its edge. 



7. Discussion 

In this work we have introduced a new, physically motivated recipe for the lateral 
expansion of the jet (in § [3]). It is based on the jump conditions for oblique shocks of 
arbitrary 4-velocity, which imply that the velocity of fluid just behind the shock front (in the 
downstream region) is in the direction of the local shock normal (i.e. perpendicular to the 
shock front at that location; /3 = n, Eq. [9]) in the upstream rest frame (which in our case 
is identifled with the rest frame of the external medium and the central source). Our new 
recipe for the lateral expansion rate of the jet (/3e ~ l/T'^Oj, Eq. [11]) has an extra factor 
of r6j in the denominator relative to the usual recipe that has been used so far {Pg ~ l/F, 
Eq. [S]). This results in slower lateral expansion relative to the usual (or old) recipe at early 
times when F > 6j, but faster lateral expansion at later times when F < 6j, i.e. once the 
lateral expansion becomes signiflcant. 

Next (in § S]) , we have implemented our new recipe as well as the old recipe in a simple 
analytic model for the jet dynamics, which is valid only for high Lorentz factors (F ^ 1) and 
narrow jet half-opening angles {6j ^1). This model shows an exponential lateral expansion 
for F < 6*0, like previous analytic models of this type. However, we demonstrate that for 
typical values of the initial jet half-opening angle (0.05 Oq < 0.2) this model is valid only 
over a very limited dynamical range for T < Oq, so that the asymptotic exponential lateral 
expansion regime is hardly reached before the model breaks down. This leads to a reasonable 



agree ment with numerical simulations over this limited range (as shown by IWygoda et al. 



201l|and in §|6]). 



This motivated us (in § E]) to generalize our relativistic model so that it would be valid 
for any values of F and 6j. This was done by switching to the 4-velocity u (instead of F) 
as the dynamical variable that we evolve (so that it would vary signiflcantly in both the 
relativistic and the Newtonian regimes), and systematically not relying on any relativistic or 
small angle approximations. Moreover, we have implemented two different assumptions for 
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the accumulation of the swept-up external rest mass, corresponding to a different variant of 
the model. The trumpet model makes the usual assumption that the working surface is the 
part of a sphere of radius R within a cone of half-opening angle Oj{R) . The conical model 
assumes that all the mass within a cone of half-opening angle Oj{R) was swept-up, so that 
once the flow becomes spherical the swept-up mass is equal to that originally within a sphere 
of the same radius. 

Our relativistic, trumpet and conical models all agree at early times when the jet is still 
highly relativistic, narrow and hardly expanded sideways (T > 9^^ 3> 1). At this stage the 
approximations of our relativistic model hold well and there are only very small differences in 
the swept-up mass between the trumpet and conical models. However, at later times when 
r < 9q^ the relativistic model enters a phase of rapid, exponential sideways expansion and 
it quickly breaks down, before becoming spherical. We note, however, that for a stratified 
or stellar wind like external medium {k = 2) the jet is closer to being spherical than for a 
uniform or ISM like external medium {k = 0; see bottom panel of Fig.^ when the relativistic 
model breaks down. 

For the trumpet and conical models, which are valid for any F or 6j, the phase of 
rapid, exponential sideways expansion largely disappears for typical values of Oq > 0.05. 
This occurs because the jet is no longer ultra-relativistic soon after F drops below 9^^, and 
once it becomes mildly or sub-relativistic its sound speed and therefore its rate of lateral 
expansion decrease compared to the ultra-relativistic regime. The conical model evolves 
somewhat faster than the trumpet model, since it accumulates external rest mass also from 
the sides of the jet, and thus it slows down faster. The smaller F results in turn in even 
faster lateral expansion rate and a larger 6j (at a given radius R or lab frame time t). 

We compared (in § E]) our analytic models to the results of 2D special relativistic hydro- 
dynamic simulations (from De Colle et al. 2011a| jbl). finding that they provide a reasonable 



description of the numerical results at all times. Therefore, they can be used for analytic 
calculations of the afterglow emission, and would provide more realistic results compared to 
previous analytic models. The main factor that significantly improves the agreement with 
simulations compared to previous analytic models is the fact that we have generalized the 
model to be valid also at modest Lorentz factors F and large jet half-opening angles 6j. 
Both our analytic generalized (trumpet and conical) models and the numerical simulations 
show that the jet first becomes sub-relativistic and only then gradually approaches spherical 
symmetry over a long time. 

For typical initial half-opening angles (^o ^ 0.05) the phase of rapid exponential lateral 
spreading is largely eliminated and it is replaced by a qua si -loga rithmic increase in 6j with 



radius R or lab frame time t. Ivan Eerten fc MacFadyenI (120111 ) have stressed that while 
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noticeable sideways expansion starts for T < 9^^, this initially involves only a small fraction 
of the total jet energy in its outer parts, and the central parts of the jet that carry most of its 
energy take longer to start spreading their energy to wider angles. While it is true that the 
jet does not remain uniform, the differences in the early growth of the angles 6f containing 
different fractions / of the jet energy, normalized by their initial value, is not very large - 
less than a factor of 2 in lab frame time or radius between / = 0.9 5 and / = 0.5, and tend t o 



become smaller for narrower ^o- This can be seen from Figure 4 of Ivan Eerten et al.l (j2011bl ). 
which also shows that as 6*0 is gradually decreased down to 0.05, its initial growth becomes 
steeper and it looks as if an early phase of expo n ential growth starts to develop, contrary 
to what is claimed in Ivan Eerten fc MacFadyenI (l201ll ). Therefore, we conclude that (i) 
although the uniform jet approximation used in our analytic models is obvious rather crude, 
it nonetheless provides a reasonable description of the energetically dominant part of the jet, 
and (ii) the prediction of our analytic models that an early exponential sideways expansion 
phase should exist for sufficiently small Oq is not only consistent with the existing simulations, 
but these simulations even show a hint for the development of such a phase. This should 
obviously be tested more thoroughly by simulations that reach even lower values of ^o- 



A phase of exponential lateral spreading was first f ound by iRhoadd (119991 ) and iPiran 



(I2OOOI ) using a simple analytic model. Later, iGruzinovl (120071 ) found a self-similar solution 
with a similar scaling. Our main conclusion (which is in agreement with lWygoda et al.ll201ll ) 
is that such a phase will occur in practice only for jets that are initially extremely narrow 
(with 6q <^ 0.05 or so), while for more modest values of 6q > 0.05 that are more typically 
inferred in GRB jets, such a phase effectively does not exist. This basically reconciles the 
long lasting apparent discrepancy between analytic models and numerical simulations. 



We thank Fabio De Colle for sharing the results of his numerical simulations. This 
research was supported by the ERG advanced research grant "GRBs" . 
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A. Appendix: Comparison to previous works 



We compare here our formulation for the jet lateral expansion rate, Eg. (|TT]). with earlie r 
work. This formula was first derived within the context of GRBs by lKumar fc GranotI (120031 ). 
who provided two different derivations. The first follows our line of argument and is based 
on the orthogonality of the shock front and the velocity of the fluid just behind it in the rest 
frame of the fluid ahead of the shock (Eq. [9]). The second derivation involves an analysis 
of the dynamical equations integrated over the radial profile. 



A result similar t o Egs. (fTOj) and (ITT1) was also rece ntly derived by lLyutikovl ( 1201 ll ). based 
on an earlier work by ShapiroHlOTQ ). Lyutikov ( 2011 ) has argued that it implies a negligible 
lateral expansion as long as F ^ 1 unless A^^ < l/F^, suggesting that with this model one 
obtains a slow sideways expansion, as seen in the numerical simulations. However, we note 
that the condition < l/F^ corresponds to Pe ~ 1. This requirement is too extreme since 
/3g ~ 1 would result in a quasi-spherical flow within a single dynamical time (since in that 
case (3d > f3r). As is well known (see also §|1]), the traditional recipe for lateral expansion 
(Eq. [8]), namely (30 ~ 1/F, already gives an asymptotic exponential growth of 9j with R 
(i.e. very rapid lateral expansion). 



The earlier work by IShapird (119791 ) discusses two possible approximations for the dy- 
namics of a non-spherical relativistic blast wave, both based on a thin shell approximation 
for the layer of shocked external medium that carries most of the energy, but with different 
ad ditional assumptions: (i) the quasi-radial approximation (used in the Newtonian regime 
by iLaumbach Sz ProbsteinI Il969l ) in which each part of the shock is assumed to move in a 
radial trajectory as if it were part of a spherical flow with the same local conditions (and 
in particular the same energy per solid angle, excluding rest energy, e = dE/dQ), and (ii) 
the iKompaneetd (jl960l ) approximation, that the pressure behind the shock is uniform, i.e. 
the same at all locations behind the shock and is proportional to the average energy density 
in the region bounded by the shock front. The first approximation assumes that the energy 
per solid angle in the flow (excluding rest energy) does not change and remains equal to 
its initial value, €{t,6) = e{tQ,6). In this sense, it basically assumes no lateral expansion 
(as the jet retains its initial angular structure in e{9) indefinitely), so that this is a model 
assumption in this case rather than a result. 



The second approximation, which was originally used by IKompaneetd (Il960l ) in the 
Newtonian regime, does not appear to be very appropriate for the relativistic regime where 
the angular size of causally connected regions is ~ 1/F ^ 1, so that that the local dynamics 
of a small portion of the flow should not be affected by the average energy per unit volume 
in the whole flow, which may b e dominated by reg ions that are not in causal contact with 
it. A simple example of how the IKompaneetd ( 1l960l ) approximation violates causality in the 
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relativistic regime is that for a uniform external medium it implies that the velocity of the 
shock front is uniforn a (i.e. depend s only on the lab frame time, but not on the location 



along the shock front; IShapirdll979l ). which necessarily implies that the flow must approach 
spherical symmetry within a few dynamical timesj§ This obviously violates causality, since 
as we discussed in the introduction, a roughly uniform jet with reasonable sharp edges cannot 
expand sideways significantly as long as T9o ^ 1, from causal considerations (since its bulk is 
not in causal contact with its edges, and it does not "know" that it is not part of a spherical 
flow and should thus start expanding sideways). 



Shapiro! (Il979[ ) reaches the conclusion that the two approximations give the same result 
in the extreme relativistic limit only because he explicitly assumed that in the quasi-radial 
approximation the energy per solid angle, e = dE/dfl, is not only independent of time, but 
also independent of the location along the shock front (this can be seen from the fact that his 
energy integral is independent of 9). This assumption quickly leads to a quasi-spher ical flow 



for a s pherical external density profile, and the non-spherical solutions obtained by [Shapiro 



( 1l979l ) arise since he considered an exponential atmosphere, which is a highly non-spherical 
external density profile. The problem of interest for us, namely the dynamics of GRB jets 
during the afterglow phase, involves a non-uniform initial distribution of the energy per solid 
angle, e(to, 9), and in such a case the two approximations are not equivalent in the extreme 
relativistic limit. Therefore, we conclude that neither of these two approximations appears to 
be appropriate for studying the dynamics or degree of lateral spreading of GRB jets during 
the afterglow phase. 



^The direction of the velocity of the fluid just behind the shock, which is along the shock normal, might 
be initially non-radial, but since the shock velocity is the same everywhere and highly relativistic, it quickly 
approaches spherical symmetry, similar to the wave left by a stone thrown into water, where the velocity of 
the surface water wave is uniform and the wave front quickly forgets the shape of the stone and becomes 
circular as its radius becomes larger than that of the stone. In our case, within a few dynamical times e 
becomes essentially independent of 6, since its local value is dominated by the recently shocked material, 
where the shock Lorentz factor is uniform. 



